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Literature 


1 .  Introduction 


Efforts  to  map  an  object  or  a  planetary  surface  usually 
include  the  reconstruction  of  the  three-dimensional  shape  of 
the  object.  There  exist  various  approaches  to  this  task: 
use  cf  overlapping  images  to  detect  and  measure  geometry 
differences  as  a  function  of  the  sensor  position 
(stereo-correlation);  use  of  brightness  changes  in  one 
image  ( photo-clinometry ,  "shape-f rom-shading" )  ;  and  ranging 
from  a  known  position  to  the  unknown  object  surface 
(profiling,  altimetry). 

One  can  denote  these  reconstructions  by 
"stereo-mapping",  since  the  greek  word  "stereo"  means 
"spatial",  "solid".  For  complex  topographical  surfaces  only 
overlapping  images  have  been  used.  Other  approaches  serve 
with  other  objects:  robot-vision  may  employ 
"shape-f rom-shading" ;  smooth  ocean  surfaces  may  be  mapped 
by  altimetry. 

Detail  in  multiple  images  for  stereo-mapping  must  be 
correlated  from  one  image  to  the  other.  This  is  done  by 
visual  matching  using  human  stereoscopic  vision,  or  a 
machine  is  used  to  perform  the  image  correlation  task. 

With  conventional  camera  photography  the  stereoscopic 
viewing  process  is  well  understood  (La  Prade,  1980).  The 
camera  is  an  analogon  to  normal  vision. 

However,  there  are  cases  where  mapping  is  done  with 
other  imaging  devices  such  as  radar  (Fig.  1.1),  e.g.  in 
adverse  climates  (tropics,  arctic)  or  in  remote  regions 
(planet  Venus).  It  is  well  known  that  stereoscopic  viewing 
can  be  used  with  radar  images  as  well.  One  also  has  an 
understanding  of  geometric  error  propagation  into  measured 
3-d  coordinates.  Only  the  limitations  of  radar  stereoscopy 
are  not  researched  to  any  extent.  Such  research  is  vital 
not  only  for  mapping  topographical  heights,  but  also  is 
needed  In  the  visual  interpretation  of  images  for  thematic 
mapping . 

Therefore,  it  is  the  purpose  of  this  report  to  review 
existing  and  to  present  new  understanding  of  radar 
stereoscopy.  For  this  purpose,  a  set  of  radar  image  data 
has  been  compiled  and  evaluated  with  respect  to  stereoscopic 
viewing.  The  range  of  data,  however,  was  found  to  be  too 
restricted  to  allow  one  to  draw  firm  conclusions.  Therefore 
a  system  for  simulation  of  radar  images  was  developed  and 
synthesized  data  were  employed  in  the  stereo  viewing 
e va lua  t ion . 


The  problem  of  radar  .stereoscopy  is  introduced  in  this 
report  by  discussion  of  stereo  viewing  that  follows  the 
extensive  body  of  work  on  the  subject  (e.g.  La  Prade,  1970; 
Leberl,  1979)*  We  then  describe  the  computer  simulation 
technique  used  to  obtain  radar  images  for  evaluation.  Both 
simulated  and  real  images  from  satellites  and  aircrafts 
serve  in  the  analysis  of  stereo-viewability  of  this  data. 

Based  on  the  conclusions  of  the  evaluation  one  now  can 
predict  the  achievable  height-mapping  accuracy.  Established 
theoretical  error  propagation  models  must  be  examined  with 
respect  to  stereo-measurability  or  viewability. 

The  conclusion  is  reached  that  numerous  flight 
arrangements  can  be  used  to  obtain  stereoscopic  radar 
images,  that  stereoscopic  viewing  often  is  only  feasible 
with  considerable  strain  on  the  observer,  and  that  height 
accuracies  can  be  obtained  for  mapping  at  small  scale,  and 
for  differential  rectification. 
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Figure  1.1  :  Radar  imaging  from  orbit,  where  s  is  the 

position  vector  of  the  sensor. 
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RADAR  STEREOSCOPY 


2 . 1  Background 


Stereo  viewing  of  overlapping  images  is  a  valuable  cool 
in  photointerpretation.  It  is  also  indispensable  for  the 
identification  and  measurement  of  homologue  points  in  two 
overlapping  images  and  subsequent  reconstruction  of  the 
three-dimensional  topographic  relief.  This  may  serve  to 
create  a  model  of  terrain  topography,  e.g.  in  the  form  of 
contours,  or  to  selectively  measure  3lopes  and  relative 
height  differences. 

This  terminology  is  based  on  the  authoritative  review 
of  LaPrade  et  al.(1980).  It  differs  from  that  used  for 
example  in  the  german  photogrammetr ic  literature  (Rinner  and 
Burkhardt  1972)  where  both  the  direct  perception  of  space 
and  the  viewing  of  overlapping  images  are  denoted  by  stereo 
viewing.  Techniques  of  generating  three-dimensional  object 
space  coordinates  from  sets  of  monocular  image  measurements 
is  the  topic  of  stereology. 

For  the  observer  radar  stereo  viewing  is  hardly 
different  from  viewing  conventional  photographs,  although 
there  exists  an  entirely  different  projection  geometry  and 
mathematical  model.  The  human  observer  perceives 
displacements  of  image  points  due  to  height  differences 
(relief  displacement).  These  form  so-called  parallax 
differences  irrespective  of  whether  we  deal  with  natural 
binocular  vision  or  with  two  images  in  stereoscopic  viewing. 

Already  a  considerable  body  of  literature  has 
accumulated  on  stereo  radar,  beginning  with  LaPrade  (1963). 
Previous  work  was  reviewed  and  some  numerical  results  were 
presented  by  Leberl  (1979).  Stereo  viewability  of  radar 
Images  was  discussed  by  LaPrade  (1970),  Graham  (1975), 
Leberl  (1975,  1978,  1979)  and  by  Kaupp  et  al.  (1982). 
Computation  of  three-dimensional  object  coordinates  from 
overlapping  images  was  analysed  by  Innes  (1964),  Rosenfield 
(  1968),  Grade  et  al.  (  1970),  Konecny  (  1972),  DBA-Systems 
(1974),  Goodyear  (1974),  Derenyi  (1975)  and  Leberl  (1972, 
1975,  1978). 

The  most  commonly  discussed  radar  stereo  imaging 
arrangements  are  shown  in  Figure  2.1. 

They  have  either  both  flights  to  the  same  side  or  each 
of  the  two  flights  at  opposite  sides  of  the  object.  Other 
arrangements  have  been  described  but  have  not  materialized: 
these  include  cross-wise  flights  (Graham,  1975)  and 
different  flight  altitudes  or  single  flight  convergent 
schemes  such  as  with  tilted  antennas  in  a  real  aperture 


radar  (Leberl,  1972;  Carlson,  1973;  Qair  and  Carlson, 
1971*,  1975  ).  It  is  not  possible  to  generate  stereo  SAR  from 
a  single  straight  flight  line  Fig.  (2.1b),  a  topic 
discussed  in  detail  in  a  preceeding  study  (Leberl,  1979). 

The  current  considerations  are  separate  from  previous 
theoretical  error  analyses.  We  discuss  the  concept  of 
exaggeration  factors  for  evalution  of  image  pairs. 


Figure  2.1a:  Flight  configurations  to  obtain  overlapping 
images  with  two  flight  lines. 


Figure  2.1b:  Flight  configurations  to  obtain  overlapping 

images  with  one  flight  line  and  real  aperture 
antenna . 


2.2  Stereo  Viewing  Considerations 
( a )  Acuity 

Measurements  in  overlapping  images  should  always  be 
made  s tereoscopically .  The  minimum  observable  retinal 
disparity  for  binocular  vision  is  reported  by  LaPrade  et  al 
(1980)  at  3"«  Muenster  (1992)  reports  a  mean  value  of  6" 
under  good  lighting  conditions-  The  optimum  can  be  achieved 
with  lines  in  object  space  that  run  parallel.  On  the  other 
hand,  there  exists  also  a  maximum  viewable  stereo  disparity 
angle:  this  amounts  to  70”  (Rinner  and  Burkhardt  1972). 
Monocularly  two  objects  can  be  distinguished  if  they  create 
an  angular  disparity  in  one  eye  of  perhaps  20".  Thus 
stereoscopy  has  a  distinct  advantage:  if  one  were  to 
monocularly  measure  the  same  point  in  two  images,  a 
measuring  error  will  be  committed  that  is  several  times 
larger  than  in  a  stereoscopic  mode. 


(b)  Viewability 

The  two  partners  of  a  stereo  image  pair  must  be  very 
similar  in  image  quality  or  thematic  content  (tone,  texture, 
etc.)  so  that  they  correlate  well,  whereas  they  should  be 
sufficiently  different  in  geometry  to  present  parallaxes  for 
height  perception.  Since  radar  is  actively  illuminating  the 
target,  differences  in  geometry  due  to  different  sensor 
positions  imply  also  illumination  differences. 

From  a  geometric  point  of  view  good  radar  stereoscopy 
therefore  seems  to  conflict  with  good  viewability.  In 
aerial  photo  interpretation  the  required  parallaxes  are 
obtained  without  any  illumination  differences  in  the  two 
stereopartners:  the  sun  illumination  hardly  changes  from 
one  photograph  to  the  next.  Stereo  viewability  is  not  a 
problem  with  photography.  It  is  the  essential  problem  with 
radar.  Figures  2.2  through  2.8  present  some  examples  of 
stereo  radar  models  from: 


(i)  aircraft  with  opposite-side  illumination  (Figure  2.2); 

(ii)  aircraft,  with  same-side  illumination  (Figure?  2 . 3 , 2 . 4  )  ; 

(iii)  same-side  illumination  from  an  aircraft  and 
satellite  SAR  (Figures  2.4,  2  .  b  )  ; 

(iv)  satellite  (SEASAT)  with  same-side  and 
opposite-side  illumination  (Figures  2.6,  2.7); 

(v)  lunar  Apollo  17  radar  with  same-side  illumination 
(Figure  2.8). 


b  km 


J 


Figure  2,2:  Aircraft  radar  pair,  opposite  side  stereo, 

of  Goodyear  Electronic  Mapping  System  (GEMS), 
over  Estrella  M  t  n  .  j  wavelength  3  cm,  s  vn  t  li  e  t  i 
aperture . 
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Figure  2.5:  Same-side  geometry  radar  stereo,  SEASAT-SAR 
wavelength  25  cm,  Flight  altitude  810  km, 
look  angles  about  ?  ^  r  . 
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Figure  2.6:  SEASAT-SAR  stereo  radar  image  pair,  flight 
altitude  810  km,  look  angles  about  20°. 
Same-side  geometry,  Los  Angeles  area. 
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Figure  2.7:  3EASAT-SAR  stereo  radar  image  pair  o  I  l.os 
Angeles,  o p po s i t o - s i d e  goometrv,  combining 
descending  orbit  with  an  ascending  one. 
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Figure  2.8:  Apollo  17  Lunar  Sounder  Experiment  (ALSE) 

image  pair  over  Montes  Appenninus  on  the  Moon, 
flight  altitude  is  112  km,  stereo  base  3.5  km, 
same  side  geometry,  wavelength  is  2  m. 

A  subjective  evaluation  of  stereo  viewability  by  an 
observer  reveals  that  influencing  factors  on  stereo  are: 

3tereo-arrangment; 
look  angles  off-nadir; 
stereo  intersection  angles; 
ruggedness  of  the  imaged  area. 

Viewability  is  thus  ensured  at  shallow  look  angles  for 
same-side  arrangements.  Opposite-side  stereo  may  be 
feasible  with  flat  or  gently  rolling  terrain.  The  limits  of 
the  actual  performance  cannot  be  defined  precisely  using 
imagery  available  today.  One  has  to  investigate  this  with 
the  help  of  an  even  larger  set  of  images,  in  particular  with 
a  larger  variety  of  cases;  image  simulation  offers  a  means 
to  evaluate  the  subjective  capability  of  an  observer  viewing 
radar  stereo  data. 
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LaPrade  (1975)  reports  on  one  experiment  with  operators 
studying  same-aide  stereo  of  flat  areas  with  man-made 
objects.  Optimum  results  were  reported  to  require  look 
angles  of  37  to  67  degrees  off-nadir  and  intersection  angles 
of  about  12  to  15  degrees.  These  intersection  angles  may 
seem  poor,  but  it  will  be  shown  later  that  radar  has  the 
potential  to  still  produce  vertical  exaggeration  approaching 
that  of  standard  photo-interpretation. 


2.3  Radar  Stereoscopic  Computations 

General  formulations  for  radar  stereo  computations  have 
been  proposed  by  Gracie  et  al.  (1970)  and  others.  The 
literature  was  reviewed  by  Leberl  (1979).  Simplified 
formulations  for  parallel  flight  lines  are  more  commonly 
employed.  For  these  recti-linear  flights  at  constant 
altitude  are  assumed  with  the  flight  direction  parallel  to 
the  object  x-coordinate  axis  (Figure  2.9). 


Figure  2.9:  Definition  of  entitles  for  radar 
stereo  computations. 


We  read  from  the  figure  that  the  object  xp  ,  yp,  2p 
coordinates  of  a  point  P  are: 

xp  =  xg  (2.1) 

y p  =  (r'2  _  r«2  +  b2)/(2B) 

zp  =  H  -  ((r'2  _  y  p  2 ) 1/2  +  (r„2  . ( B-yp  > 2 )  1  /2  ) / 2 

where  B  is  the  stereo  base,  H  is  the  flying  height. 

A  slightly  different  approach  to  compute  the  height  h 
above  a  reference  datum  is  still  with  projection  circles: 

y  =  tan  O'  (H-h) 

y  =  +  tan  0"  .  (H-h)  +  B  (2.2) 

h  =  H  -  B/  (tan  O'  +  tan  0") 

To  relate  the  object  height  h  above  a  reference  datum 
to  parallax  difference  dp  measured  on  an  image  pair  for 
given  look  and  intersection  angles,  we  have  from  Figure  9 
for  ground  range  images: 


pg' 

=  H 

cot 

V 

-  H 

cot 

0' 

pg" 

=  H 

cot 

V' 

-  H 

cot 

0" 

and 

cot  0o'  =  cot  O’  ( H-h )  /H 
cot  0Q"  =  cot  0"  ( H-h )  / H 
combining : 

dpg  =  pg"  -  pg' 

dpg/h  =  cot  0'  -  cot  8"  (2.3) 

Note  that  a  given  parallax  dpg  generates  a  different  height 
h  depending  on  Oo'and  So".  This  means  that  apparent  heights 
will  change  across  a  stereo  model,  in  contrast  to 
photographic  stereo  where  a  given  parallax  always 
corresponds  to  the  same  height  independent  of  where  in  the 


stereo  model  it  has  been  measured. 


For  slant  range  images  parallax  can  be  defined  as 
follows : 

dps  =  r"  -  p ' 

=  (H-h)  (sec  0"  -  sec  0*)  (2.4) 

Here,  unlike  ground  range  pairs,  we  find  that  zero 
heights  3till  generate  non-zero  parallaxes.  This  means  that 
the  datum  surface  in  a  slant  range  pair  will  appear  curved, 
which  we  may  calculate  by  setting  h=0,  O'  r  ©o'  and  0"  = 
Oo"  : 

dpdatum  =  H(sec  Oo”  -  sec  Oo') 

Topographic  relief  will  appear  to  lie  on  top  of  this  curved 
surface,  according  to 

dps  -  dpdatum  =  (H-h)  (sec  0"  -  sec  O') 

-  H( sec  0O”  -  sec  0Q' )  (2.5) 

We  also  know  from  the  geometry  that 

tan  O'  =  y/(H-h),  tan  0"  =  (y-B)/(H-h), 

tan  0Q '  =  y/H,  tan  0o»  =  (y-B)/H. 

These  combine  to  yield 


sec  O'  =  (1 

+  tan2 

O' 

.H2 

/(H-h)2) 

1/2 

sec  0"  =  (  1 

+  tan2 

0" 

.H2 

/(H-h)2) 

1/2 

Substituting  in  equ.  (2.5): 

(dps  -  dpdatum)/h  = 

=  ((H-h)2  +  h2  tan2  9n)1/2  _  ( ( H_h )  2  +  H2  tan2  0o')1/2  _ 
-  H.(sec  0O"  -  sec  0Q')/  h  (2.6) 


I  8 


For  H  »  h  this  simplifies  to 

(dps  -  dpdatum)/h  -  ((sec  0o?"  -  2  h/H)1^" 

(see2  0Q '  -  2h/H  )  1 _ 

sec  0o"  +  sec  0Q'  )H/h  ( ? .1 ) 


Equations  (2.3)  and  (2.7)  allow  one  to  compute  parallax 
differences  for  each  stereo  configuration  and  object  height 
difference,  and  will  be  used  for  the  evaluation  of 
exaggeration  factors  in  the  following  section. 


2.4  Stereoscopic  Exaggeration 

(a)  Definition  of  an  Exaggeration  Factor 

Stereo  imaging  geometry  is  usually  being  judged  by  two 
numbers:  The  first,  1/t,  indicates  a  ratio  between  a 
parallactic  distance  obtained  by  natural  binocular  vision 
and  the  parallactic  distance  observed  in  the  stereoscopic 
image  pair;  t  is  called  the  "total  plastic".  It  is  assumed 
that  natural  vision  and  stereo  image  generation  are  from  the 
same  distance  to  the  object.  Typically,  t  can  be  as  large 
as  several  hundred  thousand.  Since  the  figure  compares 
natural  binocular  vision  with  stereoscopic  images  it  is  not 
of  interest  in  the  current  context. 

The  second  number  is  of  greater  significance:  it 
describes  the  "flatness"  of  the  observed  stereoscopic 
impression  in  the  form  of  an  "exaggeration  factor",  q. 
Since  we  observe  stereo  imagery  with  a  central  perspective 
geometry,  we  need  first  to  understand  this  factor  for  that 
arrangement,  then  extend  the  concept  to  apply  to  radar  image 
geometry  as  well.  Stereoscopic  exaggeration,  or  affine 
stretch,  of  the  observed  object  is  a  common  concept  that  is 
discussed  in  photogrammetric  manuals  (LaPrade  et  al.  1980; 
Rinner  and  Burkhardt,  1972).  For  ease  of  reference  we 
present  this  concept. 
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Definitions  for  the  vertical  Figure  2.10b:  Definitions  for  the  vertical 

exaggeration  factor  after  La  exaggeration  factor  after  La 

Paradeetal.  (1930).  Paradeetal.  (19S0). 

Objectspace.  Stereoscopic  reconstruction. 


Figure  2.10  shows  the  central  perspective  image 
collection  geometry  for  a  pair  of  cameras  and  an  observer 
looking  at  a  stereoscopic  image  pair  through  the  lenses  of  a 
stereoscope.  The  object  is  a  pyramid  of  the  height,  hn, 
and  base  wn*  For  the  exaggeration  factor  one  defines  the 
ratio  hn/wn  the  Pyramid  in  object  space.  This 

appears  from  the  stereo  observation  as  hs/ws.  This  allows 
one  to  define  q,  the  exaggeration  factor,  as  a  measure  of 
the  flatness  of  the  observed  stereo  model: 

q  =  (hs/ws)  /  (hn/wn)  (2.8) 

It  can  be  shown  that  the  exaggeration  factor  is 
independent  of  the  dimensions  of  the  object  and  merely 
depends  on  the  so-called  "base-to-height-ratios" ,  Hn/Bn, 
in  object  space,  and  H_/B_  for  stereoviewing: 

q  =  <HS/BS)  /  (Bn/Hn)  (2.9) 

The  transformation  of  equ.  (2.8)  into  equ.  (2.9)  is 

discussed  by  LaPrade  et  al.  (1980). 

In  natural  binocular  vision  the  eyes  form  a  convergence 


angle  , 

of  about  1 4  degrees  or  less. 

This  converts  to  a 

value  for 

Hs/Bs 

of  about  4 .  One 

can  verify 

by 

experiment 
with : 

that 

optimum  stereoscopic 

viewing  is 

achieved 

Hs/Bs 

<=*  5 

so  that: 

q  «  5 

<VHn> 

(2.10) 

This  stereo  exaggeration  factor  will  now  be  applied  to 
radar  images. 

(b)  Radar  Stereoscopic  Exaggeration 

The  exaggeration  factor,  q,  as  defined  for  camera 
photography  relates  a  subjectively  observed  pyramid  in  the 
stereoscopic  model  to  the  same  pyramid  in  object  space. 
Since  we  can  relate  the  radar  stereoparallax  dp  to  the 
equivalent  photographic  stereo  case  it  is  possible  to 
compare  the  quality  of  radar  and  camera  stereoscopy.  We 
merely  need  to  find  the  photographic  ba s e- to-he igh t  ratio, 
Bn/Hn,  of  a  fictitious  camera  arrangement  that  would  produce 
the  same  parallax  dp  obtained  from  radar  for  a  given  object 
height  h.  The  exaggeration  factor,  q,  is  formed  from  equ. 
(2.9)  : 


In  case  of  a  camera  we  have  Bn  /Hn  from: 

Bn/Hn  =  dpn/hn  (2.11) 

Theref ore 
q  =  5  .  dpn/hn 

The  ratio  dp/h  needs  to  be  found  for  radar.  Using  Equ.  (2.3) 
we  find  for  ground  range  images: 

qg  «  5  .  (oot  8"  +  cot  9' )  (2.12) 

For  slant  range  images  we  may  use  equ.  (2.7)  for  the  excess 
parallax  of  topography  above  the  datum  plane: 

(2.13) 

qs  «  ((sec2  Oo"  -  2h/H)1/2  -  (sec2  Go’  -  2h/H)1/2  - 
-  sec  0Q"  +  sec  ©0*  ).5.H/h 
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V  configurations  for  the  Space  Shuttl 

fmaginp.  Radar  A  (SIR-A)  data  takes. 


We  see  in  table  2.1  that  radar  stereo  parallaxes  and 
exaggeration  factors  compare  well  with  values  obtained  from 
cameras:  as  look  angles  become  steeper,  one  has  a  more 

accentuated  stereo-effect  in  spite  of  small  stereo 
intersection  angles.  This  assumes  extreme  values  for  a  case 
such  as  Apollo  17  -  ALSE,  where  very  small  intersection 

angles  create  parallaxes  that  are  multiples  of  the  object 
height.  In  camera  photogramme  try ,  the  largest  parallaxes 
are  of  the  order  of  an  observed  height  difference,  as 
q-values  range  between  3  and  5. 

Factor  q  results  from  a  deterministic  model  of  radar 
stereo  parallaxes  and  addresses  the  question  of  an  affine 
stretch  of  the  stereo  model.  Large  exaggeration  factors 
result  from  large  parallaxes  but  do  not  ensure  high 
accuracy:  the  discussion  excludes  error  propagation  into 

measurements  of  base  width,  wn,  and  of  parallaxes,  dp. 
Theoretical  error  considerations  were  presented  by  Leberl 
(1979),  these  and  practical  results  will  be  discussed  in 
chapter  6. 


2.5  Radar  Stereoscopy  with  Non-Parallel  Flight  Lines 

Spaceborne  platforms  will  generally  not  result  in 
parallel  flight  lines  for  overlapping  images.  Instead  there 
is  convergence  of  the  flight  lines.  A  case  in  example  is 
the  Space  Shuttle  Imaging  Radar  A  (SIR  A),  as  seen  in  Figure 
2.11. 

As  a  result  one  obtains  overlapping  images  with  a 
particular  stereo  geometry  (Figure  2.12).  Stereoscopic 
viewing  is  possible;  however,  the  geometry  of  the  stereo 
case  differs  from  the  one  discussed  in  chapter  2.3* 

An  example  of  such  a  stereo  image  pair  is  presented  in 
Figure  2.13.  This  will  be  discussed  in  later  chapters  to  a 
greater  extent.  The  following  is  an  analysis  of  parallax 
geometry . 

The  two  convergent  stereo  partners  are  observed 
stereoscopically  by  placing  them  in  the  same  way  as  they 
were  generated,  i.e.  the  flight  lines  need  to  intersect 
under  an  angle, y« 

As  seen  in  Figure  2.14  one  has  to  differentiate  between 
ground  and  slant  range  presentations.  Each  image  contains 
relief  displacement,  p' ,  p".  This  combines  to  form 
observable  stereo  parallaxe  differences. 
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nvergent  ster 


Example  oi  s t e re os  cop i 
from  the  SIR -A  mission 
islands  Cephn  1  on  i  a  and 


v  a  (1  a  r 

i  in  a  y  i- 

p.i 

i  v 

l  a  k  e  n 

t  h  a  k  a  . 

•a  v  t*  r 

(  h  i 

« ■  r  * 

‘ !  r 

!•'  i  g  n  r  e  .  I  J  : 


(a)  Ground  range  presentation 


Figure  2.14  presents  the  relief  displacements  in  the 
two  images,  pg',  pg"  and  defines  a  total  parallax 
d i f f e r enc e , Apg t .  One  needs  to  attach  a  coordinate  system  to 
the  stereo  model  by  defining  an  x-axis  parallel  to  the  eye 
basis  and  y  at  a  right  angle.  The  eye  basis  is  normal  to 
the  line  of  symmetry,  ss.  Parallax  d i f f e r e nces Apg t  have 
components  in  x  and  y.  We  find: 

pg'  =  ((r'2  -(H-h)2)1/2  -  (r'2  -H2)1/2)«h.H  /  (r'2  -H2)1/2 

Pg"  =  ( (r"2  -(H-h2)2)1/2  -  (r"2  -H2 ) 1 /2 ) « h . H  /  (r"2  -H2)1/2 


Apg t  =  (pg'2  +  pg"2  -  2pg '  pg"  cosY  )1/2  2.14 
Apgx  =  (pg"  +  pg* )  siny  /  2  2.15 
Apgy  =  (pg"  -  pg’ )  cosy  /  2  2.16 

For  ground  points  not  on  line  ss  one  finds  a 
y-parallax,  pgy.  This  could  be  eliminated  by  rotating  the 
eye  basis  for  an  angle  a  : 

a  =  arctan  Apgy  /  Apgx  2.17 


Equations  2.14  through  2.17  must  be  qualitatively 
evaluated.  Fig.  2.15  presents  curves  showing  the  amount  of 
x  and  y  parallax  per  1  km  height  difference  for  various 
values  of  the  convergence  angle  y  and  look  angles  (2 '  ,  (2". 

In  the  case  of  the  Cephalonia  stereo  model,  we  find 
extreme  values  of  Apgy  or  a  as  follows: 

Apgy  max  =  I  50  m 

amax  =  18° 


These  values  can  be  neglected  for  stereoscopic  viewing. 
From  the  above  equations  of  Apgx  we  obtain  the  terrain 
height ,  h : 


h  a*  Apgx  .  tanfi  /  (2.sinY  /  2) 


2.18 


(b)  Slant  range  presentation 

As  with  parallel  flight  lines  we  find  also  in 
convergent  cases  parallax  differences  for  a  flat  area.  One 
therefore  needs  to  separate  the  observed  parallaxes  into  one 
due  to  the  actual  relief,  Aps,  and  into  the  systematic 
component,  Apo. 


Figures 

2.14  show 

that : 

po'  =  ro* 

-  (  r o'  2  _ 

h2)1/2 

=  H 

/ 

cos  n  o' 

-  H 

tan  (i  o’ 

po"  =  ro" 

-  (ro"2  _ 

h2)1/2 

=  H 

/ 

cos  0  o" 

-  H 

tan  SI  o" 

Apo 

=  ( po' 2 

+  p  o  "  2 

-  2po ' 

po" 

cosY  )1'2 

ps  • 

=  r*  - 

(  r  *  2  - 

(H-h)2 

)1/2 

=  ( H-h ) /cos  ft  ' 

-  (H-h) 

tan  ^ 

ps  " 

=  r"  - 

(r"2  _ 

(H-h)2 

)1/2 

=  (  H-h  )  /cos  ft  " 

-  (H-h) 

tan  ft 

Aps  =  (ps'2  +  ps«2  -  2ps'  ps"  cosy  J1^2 


Apsh  =  Apo  -  Aps 

We  find  a  situation  that  is  analogous  to  that  with 
ground  range  presentation.  In  the  vicinity  of  the  line  of 
symmetry,  ss,  we  can  simplify: 

h  a.  Apsh  .  cos  ft  /{  2  (  1  -  cosy  )1/2  (1  _  sin  ^  )}  2.21 

where  Apsh  results  from  EQ.  2.20  and  represents  the 
parallax  difference  above  the  deformed  reference  surface. 

In  light  of  satellite  radar  data  to  be  taken  in  polar 
regions  (antartica)  or  on  other  planets  (Venus  Radar 
Mapper),  these  considerations  promise  to  be  of  some 
relevance . 

Experimental  analysis  of  stereo  viewabllity  of  this 
imaging  configuration  will  therefore  be  included  in 
subsequent  chapters. 


3.  Image  Simulation  Approaches 


Various  algorithms  can  be  used  for  radar  imaqe 
simulation.  These  need  to  be  compared  to  develop  a 
technique  for  use  of  simulated  radar  imaqes  in  place  of  real 
ones  for  a  stereo  parameter  study.  Emphasis  therefore  is  on 
the  correct  image  geometry,  and  on  the  use  of  those 
parameters  that  are  critical  to  stereo  viewing.  Radar  image 
simulation  has  been  previously  studied  by  several  authors. 
Holtzman  et  al.  (1978)  concentrate  on  simulating 
backscatter  effects  but  neglect  geometric  effects  due  to 
topographic  relief;  Kaupp  et  al.  (1982)  concentrate  on 
geometry  for  stereo  studies,  but  neglect  the  details  of 
squint  and  flight  arrangements  different  to  the  classical 
same  side  or  opposite  side  stereo.  The  current  system 
should  allow  one  to  study  stereo  viewing  and  varying  all 
parameters  that  have  an  effect  including  squint  imaqinq  and 
different  flight  arrangements. 

The  task  of  image  simulation  may  be  split  up  into  two 
separate  tasks.  Unruh  and  Mikhail  (1981)  describe  this  as 
follows: 


"The  imaging  model  is  actually  two  separate 
models.  One,  a  geometric  model,  relates  physical 
points  in  object  space  to  physical  points  on  the  image. 
The  other,  a  spectrometr ic  model,  relates  obiect 
attributes  at  an  object  point  to  image  brightness  at 
the  corresponding  image  point." 

These  two  models  —  the  geometric  model  vs.  the 
radiometric  model,  in  our  terminologv  —  present 
transformations  between  object  and  image  spaces. 

In  the  process  of  generating  a  synthetic  imaqe  one  has 
to  relate  to  each  object  point  a  corresponding  imaqe  point 
(geometric  model) ;  two  types  of  algorithms  can  he  emploved, 
the  so-called  "Object-"  and  "Imaqe  Space  Algorithms". 
Approaches  are  discussed  in  the  following. 


3.1  Object  Space  Algorithm 


(a)  General  Description 


Simulation  is  based  on  the  Cartesian  coordinates  of  the 


points  in  a  DEM 
positions  expressed 
imaging  equations 


(Diqital  Elevation  Model)  and  on  sensor 
in  the  same  system.  After  anplving  the 
to  each  object  point  one  obtains 
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The  following  is  a  simplified  version  of  an  ob-iect 
space  algorithm.  The  input  consists  of  a  DEM  with  eauallv 
spaced  grid  points.  The  flight  Dath  is  a  straiqth  line 
perpendicular  to  the  profiles  of  the  DEM  at  a  constant 
altitude.  The  radar  beam  is  directed  at  an  ancle 
perpendicular  to  the  flight  path  (squint  anqle  t  =  0) .  All 
distances  and  positions  are  defined  in  the  same 
three-dimensional  Cartesian  coordinate  system. 
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We  find  ranqe  r  between  a  sensor  position  and  obiect 
point  P=(px>  pv,pz)  accordinq  to  Fiqure  3.2: 

r  =  SQRT { (H  -  pz) 2  +  d2)  EQ. (3.1) 


The  DEM  point  to  be  imaged  first  determines  the  sweep 
delay  s^,  the  range  r  between  sensor  position  and  anv 
other  DEM  point  along  the  profile  determines  the  imaqe 
column  coordinate  j  by 


j  =  INT  ((r-sd)/fj)  E0.  (3.2) 

where  fj  is  the  image  scale  factor  in  range  direction. 

The  image  row  i  is  determined  by  the  sensor  position  and 
thus  by  the  DEM  profile  itself. 

This  simple  approach  presents  the  user  with  such 
typical  image  geometry  effects  as  layover,  foreshor teninq 
and  radar  shadows. 

slant  range  presentation 


«  elevation  ancle  H  Fflight  altitude 

*  incidence  angle  r  slant  range 

N  Nadir  dg  sweep  delay 

RS  -  R ' S '  foreshortening  d  Distance  between  Nadir  and  P 

p’  =  c’  lay-over  Pz  Height  of  P  above  ground 

Figure  3.2  :  Imaginq  in  a  plane  perpendicular  to  the  fliqht 
path 
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(c)  Arbitrary  Fliqht  Path  and  Squint  Anqle 

However,  a  need  may  exist  to  simulate  an  imaqe  when  a 
fliqht  track  or  actual  satellite  orbit  are  known.  One  would 
have  discrete  times  and  the  cor respondinq  sensor  positions. 
The  fliqht  path  could  be  expressed  bv  an  ellipse,  hv 
splines,  polynomials  or  by  polyqons  when  simolv  connectinq 
the  given  positions.  This  serves  to  relate  any  time  t^ 
during  the  imaging  phase  to  a  sensor  position  si.  .  where 
we  express  any  sensor  position  s  in  a  local  coordinate 
system  (x,y,z) 

s=s(t)  =(sx(t)  ,sy(t)  ,s2(t) )  . 


Imaqing  is  in  a  so-called  zero-Doppler  plane  that  is 
perpendicular  to  the  corresponding  velocity  vector,  the 
first  derivative  of  s: 


*  =  S(t)  =  (Sx(t)  ,Sv(t)  ,Sz(t)  ) 

Both  £5  and  s  are  a  function  of  time  t. 

A  surface  of  constant  doppler  shift  is  in  qeneral  a 
cone.  This  is  denoted  by  a  squint  angle  t  f  0.  with  f=0 
the  cone  degenerates  to  a  plane. 

With  ground  point  P  =  (Px,Oy,oz),  a  sensor 
position  j;  with  the  squint  anqle  t,  and  velocity  s,  we 
find  from  Fiqure  3.3  that  a  relation  exists  between  o,  s, 
*  and  t: 

- 

_  * -  =  cos  ( 90°-t)  =  sin  l  EO.  (3.3) 

|s|  |p-s| 

Assuming  s  and  js  to  derive  from  polynomials  in  t, 
this  equation  can  be  solved  for  everv  point  P  of  the  around 
to  be  imaged.  The  solutions  are  the  set  of  { t^ , 

k=l,...m}.  It  is  a  mathemetical  problem,  that  the  closer 
the  flight  path  is  approximated  by  the  polynomial,  the  more 
it  oscillates.  This  leads  to  a  quantativelv  biq  set  of 
solutions,  from  which  a  subset  of  possible  solutions  ha^e  to 
be  extracted.  This  is  no  trivial  task.  T’his  subset  f  t^f 
1=1,..., n  |  n  <  m}  determines  sensor  positions 

From  a  specific  position  s  at  time  t  the  slant  ranoe  r 
i  s 

r=  i p  -  s |  FO. (3.4) 

This  is  displayed  in  a  slant  ranqe  presentation.  Tn  a 

ground  range  presentation  one  finds  the  column  coordinate  a 
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as 


g  =  SORT { r 2  -  H2)  EO. (3.S) 

where  H  is  meant  to  be  the  altitude  of  the  sensor  at  imaqinq 
time  t.  It  is  thus  the  current  value  at  the  time  of  fliqht. 

Time  ( t)  defines  the  image  row  (azimuth  coordinate), 
whereas  (r)  respectively  (q)  define  the  imaqe  column  or 
range  coordinate.  (See  EQ.  3.2) 


3.2  Image  Space  Algorithm 


(a)  General  Description 

The  image  space  approach  starts  from  equidistant  imaqe 
coordinates  (i,j)  of  the  output  imaqe.  Radar  imaqinq  time  t 
and  range  r  can  be  derived  from  these  coordinates  and  define 
a  projection  circle.  The  cor respondinq  object  point(s)  in 
the  DEM  must  be  found  bv  intersectinq  the  circle  with  the 
DEM.  The  intersectinq  points  are  usually  not  exact  DEM  qrid 
points;  the  intersection  of  arc  and  ground  is  thus  found  hv 
interpolation.  This  kind  of  approach  is  useful  when  object 
space  points  must  be  associated  with  previously  specified 
image  locations  directly  (See  Greeve  and  Goonev,  1974) . 


(b)  Deriving  sensor  position  s  and  ranqe  r  from  imaqe 
coordinates 

The  coordinates  (i,j)  in  the  output  image  serve  as 
input  values  for  the  simulation  oroaram,  where  (i)  is  a 
function  of  time  (t)  and  (j)  is  a  function  of  slant  ranoe 
(r)  or  ground  range  (g) .  We  find: 


t  =  (  i-1)  *f .  +  tQ 
r  =  ( j-1) *f j  +  sd 
g  =  (j-1)  *f .  +  s<a 
r  =  SORT (H2  +  g2) 


where 


EO. (3.6a) 
EO. (3.6b) 


EO. (3.6c) 


t.... imaging  time  of  coordinate  (i,j) 
fc0... start  time  of  imaging 
£^... scale  factor  in  azimuth-direction 
fj... scale  factor  in  ranqe-direction 


The  platform  position  S  (described  bv  vector  s)  and  the 
corresponding  velocity  vector  x  can  be  calculated  as 
functions  of  t. 
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In  addition  one  has  an  antenna  (sensor)  coordinate 
system  (u,v,w)  with  the  oriqin  at  the  platform  position  an^ 
u  along  the  platform  velocity  vector.  P  has  the  coordinates 

£  =  ( Pu r Pv r  Pw) t  the  relationship  between  platform 

position  and  point  P  to  be  imaqed  in  the  local  coordinate 
system  is  defined  by  EQ.  (3.3). 

The  antenna  coordinate  system  is  defined  bv  three  unit 
vectors  jj,v,w: 

u  =  s/| s| 

v=(sTS)/s*S  EQ.  (3.7) 

w  =  ( u  *  v)  /  £  *  v 

EQ.  3.3  degenerates  to 


u  *  p  =  r  *  sint  eq  .  (3.3) 

in  this  coordinate  system. 

Introducing  the  elevation  angle  0  the  followinq 
relationships  exist  between  the  coordinates  P^ f  p 
Pw,  elevation  angle  0  and  the  squint  anqle  t: 

Pu  =  r  *  sin  t 

Pv  =  r  *  SORT ( s in  w  -  sin^t) 

Pw  =  -  r  *  cos0  EQ.  (3.9) 

EQ. (3.3),  EQ. (3.9)  and  EO. (3.8)  are  alternative 
formulations  of  the  relationship  between  the  existina 
parameters. 

As  both  of  the  coordinate  svstems  are  cartesian,  the 
conversion  between  the  local  and  the  antenna  svstem  consists 
of  a  shift  and  rotation  with  shift  vector  q  and  a  rotation 
matrix  R. 

(pu'Pv'Pw}  =  R*( (PX'Pv'pz>  “  2  )  iQ .  (3.10a) 

or  in  the  other  direction: 

(PX'DyrPz)  =  HT* (Du,Pv,pw)  +  q  EQ-  (3.10b) 

(  q  is  of  course  found  to  be  s  and  R  to  be  (uT, v^, wT) , 
alT  expressed  in  the  local  Cartesian  svstem). 

is  the  transposed  matrix  of  R  . 


(d)  Ranqe  Sphere  and  Doppler  Cone;  Intersection  Circle 


Any  imaqe  point  creates  thus  the  followinq  qeometrv; 

Time  t  (derived  by  the  imaqe  row  coordinate  i)  defines 
the  sensor  location  S  and  its  velocity  vector  $5. 
(s,S;t)  determine  a  cone  (Doppler  cone) ,  whose  peak  is  in 
S,  and  all  generating  lines  form  an  angle  of  (90  -f)  with  s. 

The  range  column  j  defines  the  slant  range  r  and 
creates  a  sphere  with  the  radius  r  around  its  center  S. 
Intersection  of  the  sphere  and  cone  define  a  circle  as 
illustrated  in  Figures  3.5  and  3.6.  We  define  the  sphere 
(s, r)  ; 


u2+v2+w2  =  r2 
The  cone  (s;u;t)  is 


EQ.  (3.11) 


(v2+w2) /r2cos2f  -  u2/r2sin2t  =  0  EQ.  (3.12) 

The  intersection  of  cone  and  sphere  defines  a  circle: 

v2  +  w2  _  r2*cos2t  EQ.  (3.13) 

in  the  plane 

u  a  r  *  sinf. 

Each  image  point  then  represents  a  circle  that  needs  to 
be  intersected  with  the  qround.  In  the  simulation  task, 
this  is  an  intersection  between  circle  and  DEM. 
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Figure  3.5:  Doppler  cone  and  range  sphere. 
Perspective  view. 
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Figure  3.6:  Doppler  cone  and  range  sphere. 
Vertical  view. 


(e)  Intersectinq  the  Projection  Circle  with  the  DEM 


The  intersection  between  circle  and  DEM  is  found  in  two 
steps.  Let  a  ray  from  S  form  an  anqle  9  with  the  nadirline. 
Its  intersection  with  the  projection  circle  produces 
position  £(0)  =  (Pu»Pv»  Pw)©  using  EO.  (3.7).  In  a 
next  step  p  is  converted  to  £  =  (Px,Py»Pz)0  usinq 

EQ.  (3.9d) .  Heiqht  pz)  can  be  comoared  with  real  heiqht 
z  at  planimetr ic  position  (Px,pv)  where  z  is 

approximated  linearly  by  the  fourvDEM  corner  points.  Bv 
changing  0  monotonously  increasinq  or  decreasing,  the 
intersection  points  will  be  below  or  above  the  actual 
terrain;  the  canqe-over  from  below  to  above  or  vice  versa 
indicates  an  intersection  within  an  interval  [0.  p. 

i » ” i— li • 

Then  this  new  interval  is  subdivided  and  the  computations 
are  continued  to  define  the  intersection  point  more 
precisely.  Finally  the  DEM  segment  [0^,0.  is 

approximated  by  a  plane  and  intersected  with’ thi  circle. 

Each  circle  has  to  intersect  the  DEM  at  least  once,  but 
more  than  one  solution  is  possible  and  meaninqful.  T’he 
occurrence  of  more  than  one  solution  is  due  to  the 
"lay-over";  more  than  one  DEM  point  is  to  be  imaqed.  In 
other  words:  more  than  one  ground  point  has  a  qiven  ranee  r 
from  a  platform  position. 


(f)  Arbitrary  Flight  Path  and  Squint  Anqle 


The  approach  used  for  image  simulation  with  the  imaqe 
space  algorithm  defines  a  flight  path  in  one  of  several 
forms:  lookup  tables,  polynomials,  splines  etc.  It  is 
mereley  the  problem  of  relating  the  imaqe  time  to  the 
corresponding  sensor  position.  This  allows  one  to  enter 
flight  perturbations  (pitch,  yaw)  and  to  study  their  effect 
on  the  image. 

Imaging  with  t  =  0  as  well  as  f  =  ft  is  possibly. 
However,  for  the  latter  case  there  is  a  need  to  additionally 
rectify  the  image. 


(g)  Geometric  Rectification  of  Squint  Angle  Imaqes 


All  points  imaged  from  one  sensor  position  are  found  in 
one  image  line.  In  the  case  of  t  =  0  the  corresponding 
object  points  were  lined  up  in  a  plane  perpendicular  to  the 
flight  path  and  distortions  appear  in  the  ranqe  direction 
resulting  from  the  height  variation  on  the  ground.  Tn  the 
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other  case,  however,  the  imaqed  points  appear  in  the  obiect 
space  as  elements  of  a  cone  (the  Doppler  cone)  and  thus  the 
image  pixel  locations  have  to  be  rearranged  as  to  where  thev 
would  appear  in  the  case  of  imaging  with  t  =  0.  The 
rectification  is  explained  in  Fig.  3.8  for  a  forward 
looking  squint:  All  intersection  points  of  DEM  and  cone 
appear  in  one  image  line,  as  they  are  imaqed  from  the  same 
sensor  position  S. 


Figure  3.7:  Rectifying  "squint  angle"  imaqes  qeometr icallv: 
new  sensor  position  S^'  and  new  ranoe  (r*) 
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All  intersection  points  of  DEM  and  cone  appear  in  one 
image  line,  as  they  are  imaqed  from  the  same  sensor  position 

Si  (Si  =  peak  of  cone  at  imaginq  time  ti,  where  ti 
is  the  time  of  imaging  row  i) .  The  location  where  the  image 
of  one  DEM  point  P  lies  within  this  line  is  dependent  on  the 
distance  r  between  and  P.  Without  squint  anqle  point  P 
would  have  been  imaged  from  another  position,  r  *  sin  t 
different.  Considering  the  scale  factors  in  azimuth  and 
range  direction,  imaqe  point  P'  should  he  moved  forward  to 
this  new  image  line.  For  a  slant  ranqe  presentation  the  new 
range  r*  =  r  *  cos  t  determines  the  imaqe  column  of  the 
rectified  image. 

Fig.  3.8  show  a  simulation  result  of  usinq  a  squint 
angle  £  =  30°.  Other  parameters  for  this  simulation  are 

Flight  altitude  =  264  km 
Elevation  anqle  to  near  range  =  70°. 


Text  for  Figures  3.8  (a),(h)  and  3.9  on  next  paqe: 

Fig.  3.8  (a)  is  the  raw  simulation  result.  Fig.  3.8 
(b)  shows  the  imaqe  after  rearranging  the  imaqe  pixel 
locations  to  the  geometry  of  t  =  0.  The  radiometric 

distortions  caused  bv  the  squint  cannot  he  rectified.  For 
comparison,  Fig  3.9  is  a  presentation  of  a  simulation 
without  squint  angle,  hut  with  exactly  the  same  other 
imaging  parameters. 
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3.3  Radiometric  Model 


Once  the  relationship  between  imaqe  and  obiect 
coordinates  is  established  one  can  proceed  to  assiqn  density 
numbers  to  each  location.  These  values  will  vary  between 
black  (shadow)  and  white  (tarqet  reflectinq  specularly 
towards  the  sensor) . 


(a)  Shadow 


The  determination  of  areas  in  radar  shadow  is  a  hidden 
surface  problem.  Hidden  surface  elimination  is  well 
documented  (e.g.  Sutherland  et  al,  1974)  .  '’’he  problem 

remains  to  find  a  fast  alqorithm  for  each  special 
application. 

Hidden  surface  computations  always  implv  some  kind  of 
sorting.  In  the  case  of  an  object  space  alqorithm  movinq 
along  the  DEM  profiles  from  near  to  far  ranqe,  the  array  is 
already  presorted  and  shadow  areas  can  be  easily  found  bv 
the  following  method. 

For  each  point  P  the  incidence  angle  91  (defined  as 
the  angle  between  the  surface  normal  .at  P  and  the  line 
between  sensor  and  P)  is  calculated.  If  91  passes  90° 
in  P,  then  a  straight  line  from  the  sensor  throuqh  point  P 
is  used  to  compute  intersection  point  O.  The  area  from  P  to 
Q  is  in  the  shadow.  At  0  the  imaging  computations  continue. 

Other  algorithms  prove  to  be  more  difficult  to  use 
since  no  information  is  available  at  imaginq  time  if  a  DEM 
point  is  hidden  by  other  parts  of  the  DEM  or  not. 

Therefore  the  shadow  is  computed  in  a  post-processinq 
step.  Imaging  results  in  information  on  corresponding  imaqe 
pixels  and  DEM  cells  such  as  elevation  angles  and  distances 
from  nadir  line.  These  must  be  sorted  with  respect  to  the 
corresponding  time  rows  and  ranqe  columns.  Sortinq  in  time 
rows  is  needed  since  DEM  cells  are  only  hidden  by  other  DEM 
cells  that  are  seen  from  the  same  sensor  position. 

Sorting  along  the  range  column  is  in  two  steps.  First, 
points  are  sorted  with  increasinq  qround  distance  from 
nadir.  Second,  a  check  is  performed  on  the  elevation  angles 
for  strictly  increasinq  values.  For  non-increasinq  values 
of  the  elevation  angle  the  corresponding  imaqe  pixels 
represent  radar  shadow. 
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The  imaginq  technique  always  results  in  an  almost 
presorted  data  field.  With  "Quicksort",  (Hoare,  1962)  one 
of  the  classical  very  fast  sorting  algorithms,  the  oomputina 
time  would  be  0  (n  *  log^n)  for  a  best  case,  where  0  means 
"in  the  order  of",  and  n  is  the  number  of  elements  to  sort. 
It  is  not  applied  because  of  its  property  to  increase 
computing  time  for  presorted  arrays  (Maurer,  1974). 
"Heapsort"  (Williams,  1964)  was  applied  instead;  this 
method  has  a  worst  case  of  0  (n  *  log2n)  once  the  field  is 
changed  into  a  heap,  and  proves  to  be  negligible  in  rpij  time 
for  this  application. 


(b)  Assigning  Grey  Values 


The  grey  value  to  be  assigned  to  an  image  coordinate 
(i,j)  is  a  function  of  the  incidence  angles  Q1  appearing 
at  the  corresponding  ground  points.  The  angle  enters  into  a 
backscatter  function  <p ( © 1 )  . 

The  incidence  angle  is  defined  as  the  angle  between  two 
vectors,  namely  the  normal  onto  the  HEM  and  the  incidence 
ray,  (£-s)  (ray  from  sensor  position  S  to  qround  point  P) . 
In  the  current  context  of  stereo  studies  it  was  considered 
sufficient  to  use  one  backscatter  function  for  an  entire 
image.  One  could,  however,  define  different  such  functions 
depending  on  planimetric  positions  x,y  in  the  DEM.  mhis 
would  allow  one  to  differentiate  between  e.a.  forest,  sand, 
water,  ice,  etc.  However  this  would  imply  the  existence  of 
thematic  information  in  addition  to  the  mere  topographical 
data  in  the  DEM. 


So  far  standard  backscatter  functions  of  Hagfors  (1969) 
or  Muhleman  (Kobrick,  1982,  personal  communication)  and  a 
cosine  function  have  been  used  accordina  to  EO's.  3.14: 

1 


^HAG 


logio  sqrt  t 


(cos4©1  +  C* 


sin^Q) 


Eg.  (3.14a) 


^MUL  =  l°9l0 


cos©  i*A3 


(sin©1  +  A*cos©i) 3 


EO.  (3.14b) 


^COS 


=  cos©1 


EO.  (3.14c) 


where 


C  ...  is  Hagfors'  constant 
A  ...  is  Muhleman's  alpha 
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The  resulting  (J)-values  must  he  converted  to  diqital 
density  numbers  DN. 


3.4  Discussion  and  Choice  of  Implementation 


The  algorithm  design  was  the  result  of  experiments  with 
a  simplified  simulation  using  the  object  space  approach. 
Shadowing  can  be  included.  This  type  of  simulation  is  a 
useful  and  fast  tool  if  no  geometric  perturbations  must  be 
studied. 

A  more  generally  applicable  approach  is  with  the  use  of 
an  image  space  method.  It  allows  for  flexible  consideration 
of  perturbed  fliqht  paths,  resulting  in  shadows,  layores 
etc.  It  actually  models  the  imaging  process. 

Benefits  from  this  flexibility  accrue,  for  example  in 
the  SIR-A  case: 

SIR-A  was  not  designed  to  deliver  stereo  imaqe  pairs. 
Images  taken  over  the  same  area  are  with  the  same  look  angle 
(app.  50°  at  mid  swath)  but  with  crossinq  flight  tracks. 
This  proved  to  be  valid  stereo  imagery  (see  chapter  2) .This 
type  of  imaging  can  be  studied  with  a  simulation  method  that 
offers  the  flexibility  of  arbitrary  fliqht  directions. 

The  actually  implemented  technique  was  also  required  to 
allow  for  investigations  regarding  squint  anqle 
arrangements. 

One  therefore  has  now  available  a  method  of  simulation 
that  can  use: 

-  available  flight  recordings, 

-  arbitrary  resolution  and  size  of  DEM, 

-  arbitrary  resolution  and  size  of  simulated  image, 

-  squint  angle  imaging, 

-  range  processing  (ground  resp.  slant  ranqe  geometry) , 

-  choice  of  elevation  anqle  (or  sweeo-delav) , 

-  choice  of  look  direction  of  SLR, 

-  uniform  backscatter  curves  accordinq  to  EO.  3.14  (b) . 

In  addition  there  exist  a  simple  method  that  is  faster 
and  applicable  where  special  parameters  need  not  be  studied. 

The  product  of  the  software  implementation  is  named 
SIMRISA  (SIMulation  of  Radar  images  using  an  Imaqe  Space 
Algorithm)  and  described  in  appendix  A.  Emphasis 'was  put  on 
modularity  of  the  system.  The  proqram  can  be  easily 
changed,  e.g.  additional  thematic  input  data  sets  could 
allow  for  use  of  different  backscatter  curves  or  use  of 
different  look-up  tables  within  one  DEM  area. 


4.  Verification  of  the  Simulation  Proqram 
4.1.  Using  Geometric  Shapes 


Synthetic  DEM's  serve  to  demonstrate  the  features  of 
the  simulation  system.  Fiqures  4.1  (a),  (b)  show  one  of 

these  synthetic  models.  The  two  profile  lines  anpear  again 
in  Fig.  4.2,  where  the  slant  range  geometry  for  these 
special  cases  is  explained. 

The  illumination  (or  problem  of  hidden  surfaces)  is 
presented  in  Fig.  4.3.  Here  as  well  as  in  the  followino 
radar  images  of  this  model  the  ground  level  is  not 
illuminated  (or  imaged)  but  set  to  a  constant  value  to 
enhance  the  shadows  and  geometric  features  of  the  imacie. 

The  flight  path  is  chosen  to  show  distinct  lav-overs 
and  shadows. 

Fig.  4.4  presents  a  synthetic  slant  range  oresentation 
of  the  data. 

An  assumption  of  diffuse  scattering  is  made.  with 
varying  incidence  angles  on  a  flat  ground  from  near  ranae  to 
far  range  the  grey  values  change,  too,  from  hiqh  values  to 
low  ones.  The  geometric  shapes  change  in  both  their  imaqe 
geometry  as  well  as  their  brightness.  Verv  high  return 
values  indicate  layovers  (more  than  one  return  value  arrives 
at  the  receiver  at  the  same  time) .  Grev  value  0  indicates 
radar  shadow.  Not  onlv  the  geometry  but  also  briqhtness  of 
an  object  changes  as  it  is  imaged  at  different  look  angles. 
This  is  a  main  cause  of  difficulties  in  stereo  radar. 

Fig.  4.5  (a)  and  (b)  present  an  opposite  side  stereo 
pair.  Here  the  parameters  were  chosen  like  thev  might  be 
with  the  ERS-1,  at  a  flight  altitude  of  777  km  and  an 
incidence  angle  of  23°  at  the  center  of  the  (80  km  wide) 
swath.  However,  stereo  viewability  is  not  possible  for  this 
specific  terrain  and  sensor  configuration.  (The 

configuration  is  presented  in  Fig.  4.6) 

The  data  set  can  serve  for  educational  purposes  since 
it  allows  one  to  follow  through  a  series  of  parameters  with 
the  same  known  shapes. 


Figure  4.2:  Imaging  geometry  for  slant  range 
presentation  corresponding  to  the 
profiles  in  Figure  4.1  (b). 


Figure  4.3:  Radar  illuminated  height  model 
(shadow  and  illuminated  areas 
can  be  clearly  separated). 


Figure  4.4:  Radar  slant  ran 


g  e  presentation. 


4.2  Using  SIR-A  Data 


The  Shuttle  Imaging  Radar  A  experiment  was  carried  out 
in  November  81.  The  swath  width  was  approximately  50  km, 
with  an  incidence  angle  of  50°  at  mid-center.  The  wave 
lenqht  was  23  cm  (L-band) . 

(a)  Greek  Islands 


A  well-suited  data  set  for  radar  stereo  studies  was 
generated  during  the  SIR-A  mission  over  two  of  the  Ionian 
islands,  Cephalonia  and  Ithaka.  The  islands  are  situated 
between  38  00'  N  and  38°30'  N  latitude  and  between 
20°19'  e  and  20°49'  E  longitude. 

For  stereo  studies  a  DEM  had  to  be  generated.  This  was 
done  using  a  topographic  map  at  scale  1:200  000  in  an  U^M 
projection  with  contour-lines  at  a  200  m  interval.  The 
countour-lines  were  digitized  into  a  qrid  with  a  mesh  size 
of  1"  in  latitude  and  longitude. 

Cephalonia  is  the  largest  (743  sqkm)  island  of  the 
Ionian  group  and  also  the  most  mountainous  one.  The  highest 
elevation  -  mountain  Aenos  covered  with  firs  and  stone-pines 
-is  1628  m.  Population  is  mainly  on  the  shores,  like  the 
capital  Argostolion  (7100  inh.),  Lixuri  (3400  inh.)  or  Sami 
(1000  inh.).  The  vegetation  is  rich  and  nearly  tropical. 
The  economy  is  based  mainly  on  vineyards  and  olivetrees. 

Ithaka  -  native  island  of  Odysseus  -  is  rocky  and  not 
very  fertile.  It  is  divided  into  two  mountainous  areas 
connected  by  a  land-bridge  only  600  m  wide.  The  hiqhest 
elevation  is  mountain  Nitros  (806  m)  in  the  northern  part. 
Ithaki  is  the  main  city  (2300  inh.) 

During  the  SIR-A  flight  Cephalonia  and  Ithaca  were 
imaged  two  times:  on  Data  Take  32-33  and  Data  Take  37A. 

The  two  orbits  crossed  at  an  angle  of  approximately  34° 
and  created  two  images  sufficient  to  produce  observable 
stereo  parallax.  The  stereo  pair  was  shown  in  Fiq.  2.13. 
The  swath  width  is  50  km  ,  look  anqle  47°.  SIR-A  was 
optically  processed. 

The  radar  images  show  all  main  cities  in  Cephalonia  by 
their  bright  return  values:  Arqostolion,  Lixuri,  Sami  and 
Karavomilo  from  the  map  in  Fiq.  4.7  can  be  identified  on 
the  images.  The  mountainous  region  in  the  central  part  of 
Cephalonia  (forested  reqion)  is  a  quite  homogeneous,  diffuse 
scatterer  and  shows  the  brightness  depending  on  the 
incidence  angle.  The  runwav  of  the  airport  south-west  of 
Argostolion  is  a  specular  scatterer  and  returns  no  echo  to 
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The  simulations  of  Fiq.  4.8  are  based  on  a  homoqeneous 
backscatter  curve  for  the  whole  test  area.  The  diffuse 
scattering  of  the  forested  central  part  can  be  modeled  bv 
the  Hagfors  law.  Comparison  of  the  reai  and  simulated  imaqe 
in  the  flat  western  peninsula  shows  the  limitation  of  the 
simulation  if  no  thematic  information  is  used.  Also  the  DEM 
has  accuracy  limitations  so  that  smaller  features  are  lost 
in  the  simulation. 


s  I  «  -  0  37A 

HAC 


Figure  4.8  :  Simulations  of  Greek  Islands,  usinn  STR-A 
Parameters.  (C im ino/Elach i ,  19R?) 


(b)  Northern  California 


Another  test  site  was  found  in  North  California 
(latitude:  40°  -  41°  N,  longitude:  123°  -  121°  E) . 

The  elevation  model  was  obtained  from  the  U.S.  Geological 
Survey.  It  is  in  the  U.T.M.  projection  and  is  obtained 
from  contour-lines  on  a  1:250000  topographic  mao. 

To  achieve  a  resolution  comparable  to  the  existinq 
radar  image  the  DEM  was  reduced  -  each  4  heiqht  points  of  a 
grid  cell  were  averaged  into  one  value,  resulting  in  a  DEM 
cell  diameter  of  126  m.  Usinq  a  disolav  screen  showinq 
512x512  pixels  this  test  site  represents  a  square  area  with 
approximately  64.5  km  on  each  side. 

The  corresponding  real  radar  image  was  obtained  hv 
SIR-A,  orbit  21,  data  take  24B.  The  resolution  of  the  imaqe 
was  roughly  40  m;  originally  SIR-A  images  were  analog;  the 
test  data  were  digitized  with  a  pixel  diameter  of  about  25 
m.  To  cover  a  comparatively  larger  area  in  a  512x512  pixel 
array,  each  25  image  points  were  averaged  into  one  to  obtain 
about  the  same  pixel  resolution  as  for  the  simulated  image. 

The  topo  data  file  was  preprocessed  for  the  earth 
curvature  and  rotated  to  be  parallel  to  the  direction  of  the 
Space  Shuttle.  This  way  the  simplified  version  was  used  for 
simulation.  Ultimately  a  slant  ranqe  presentation  was 
obtained. 

Real  and  simulated  images  of  the  test  area  appear  in 
Fig.'s  4.9  (a)  and  (b) ;  a  backscatter  function  after 
Haqfors  (see  3.2)  served  in  the  simulation.  The  results 
appear  to  be  sufficiently  realistic  to  serve  in  stereo 
viewability  experiments. 


Work  for  the  simulation  of  the  North  California  test  site  was 
done  while  one  of  the  cooperators,  G.  Domik,  was  a  visiting  research 
affiliate  at  the  Jet  Propulsion  Laboratory,  Pasadena,  California.  nr. 
M.  Kobrick  was  the  supervisor  during  the  visit.  His  support  is 
gratefully  acknowledged. 


Figure  4.9  :  Real  and  simulated  image  from  N-Oalifornia 

(a)  SIR  A  :  Orbit  21  /  Data  T’ake  24  r. 

(b)  Simulation  usina  SIR-A  Parameters. 
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Just  west  of  the  presented  area  lies  the  center  of 
another  swath  crossinq,  produoinq  a  cross  anale  stereo  pair 
(crossinq  angle  14°).  Both  the  real  and  the  simulated 
image  are  presented  in  (Kobrick  et  al,  1RR3). 


4.3  Using  SAR-580  (Oetztal). 


The  SAR-580  Experiment  55-0  "Snow  and  Ice  Parameters 
from  Radar  Data"  was  associated  with  a  qlacierized  reqion  in 
the  Austrian  Alps  (Rott  H.,  1983).  T’he  main  obiectives  were 
directed  at  applications  of  SAR  data  for  mapping  of  snow  and 
glaciers  in  mountain  regions. 


On  Julv  7,  1981,  one  pass  (no. 149)  was  acquired  in  X- 
and  C-Band  (both  HH-Polar isation) ,  coverinq  a  ground  swath 
of  13  km  length.  Fig.  4.10  shows  a  sketch  map  of  the 
ground  swath  covered  by  the  digitally  processed  SAR  data; 
it  is  centered  at  about  46.8°  n  and  10.9°  p.  The 

testsite  covers  an  altitude  ranqe  between  2000  m  and  3500  m 

a .  m .  s .  1 .  , 

I 

1 


4.11  :  Digital  elevation  model  of  Oetztal.  The 
actual  radar  imaqe  is  only  Dart  of  the 
data  displayed  above. 

above:  axonometric  presentation  of  heiqht  model, 
below:  Illuminated  heiqht  model. 
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In  Fig,  4.12  the  imaqinq  parameters  of  the  profiles  in 
Fiq.  4.10  are  explained.  Simulation  was  carried  out  usinq 
these  parameters.  Also  the  corrections  from  a  slant  range 
dependence  to  a  ground  ranqe  presentation  for  a  flight 
altitude  of  6100  m  a.m.s.l  to  a  nominal  altitude  of  2600 
were  m  a.m.s.l  carried  out  during  processing.  Comparisons 
of  the  real  and  simulated  image  in  Fiq.  4.13  (a)  and  (b) 
show  the  same  geometry.  The  radiometry  is  not  matched  to 
the  real  one  since  this  would  lead  beyond  the  scope  of  the 
current  research. 


Figure  4.12:  Surface  relief  at  the  profiles  1  and  2  of  the  test 

site  0-1  with  radar  imaginq  geometrv.  Zones  of  radar 
glaciers  are  hatched.  0^  =  50°,  ©2  =  55°, 

©3  =  60°,  ©4  =  650,  o5  =  70o 
antenna  incidence  angle. 
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flight  path 


0  E  T  2  T  ft  L 
SIHUUTED  RADAS  IMAGE 


5.  Discussion  of  Stereo-Viewabilitv 


Evaluation  results  of -stereo-viewabilitv  are  discussed 
showing  tables  for  "same  side",  "same  side  with  applied 
squint  angle"  and  "cross-angle"  stereo  .  Different 
investigators  were  requested  to  rate  the  viewabilitv: 
Persons  without  previous  stereoscopic  viewinq  and  measurinq 
experiences,  and  two  photogrammetr ists  with  such  experience. 
The  viewability  is 

-  described  by  a  number  between  1  and  10  by  the  first 
group,  where  1  denotes  "no  viewabilitv"  and  10 
"excellent  viewability".  There  is  also  a  distinction 
between  1-4  and  5-10,  the  lower  numbers  indicatinq  "not 
accaptable"  and  the  higher  ones  "acceptable"  stereo. 

-  described  by  a  letter  q,  m  or  n  (by  the 
photogrammetr ists) ,  where  g  means  "good  visible",  m 
means  "marginal"  and  n  means  "not  visible". 

This  second  grading  shows  stereoscopic  impression  to  a 
much  higher  extend. 

The  model  area  for  these  tests  was  the  Greek  Island 
test  site,  which  can  be  characterized  as  a  "hiqh  relief 
area".  The  site  is  described  in  more  detail  in  4.2  (a). 
According  to  the  stereo  viewing  considerations  (chapter  2.2 
(b) ) ,  the  two  stereo  partners  should  be  similiar  in  thematic 
content.  Thus  a  homogeneous  backscatter  curve  was  used, 
homogeneous  in  the  sense  of  not  differentiating  between 
thematics  in  the  image  as  well  as  in  the  sense  of  using  the 
same  backscatter  curve  for  a  whole  set  of  stereo  imaqes  to 
be  evaluated. 


5.1. 


Same  Side  Stereo  —  Different  Elevation  Anqles 


Table  5.1  shows  the  quality  of  the  viewahilitv  of  same 
side  stereo  for  twelve  simulated  images  (all  in  slant  ranqe 
presention) .  Their  elevation  anqles,  also  denoted  as  "look 
angles",  range  from  10°  to  80°  (See  Fiq.  5.1  for  stereo 
configurations  and  Fig.  5.2  for  examples).  Stereo  v’ewinq 
for  the  northern  part  of  the  island  (heiqhts  less  than  1000 
m)  is  possible  for  a  much  wider  ranqe  of  intersection  anqles 
than  viewing  the  southern  part  of  Ceohalonia,  where  the  area 
is  very  rough  and  higher  and  thus  creates  more  shadow  for 
shallow  look  angles  and  more  foreshortening  and  layover  for 
steep  look  angles. 

The  highest  ranked  stereo  pairs  (evaluated  9-10/q)  are 
closer  investigated  (table  5.2).  An  intersection  anqle 
between  10°  and  25°  was  found  to  be  best. 

The  photogrammetr ists  found  60°  intersection  anqles 
still  to  be  viewable.  The  back-slopes,  however,  were  nudged 
to  be  more  difficult  to  fuse  than  slopes  facing  towards  the 
antenna. 


DEM  Digital  Elevation  Model 

SF  ....  Satellite  at  first  flight 

SL  ....  Satellite  at  last  flight 

a . Altitude 

s.d.  .  ...  Sweep  Delay 
. Elevation  Angle 

. 10* 

«L . 80* 


Figure  5.1:  Configurations  of  same  side  stereo 
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Same  Side  Stereo 


Table  5.1:  Same  side  stereo  evaluation. 


RANK  ORDER 

ELEVATION  ANGLES 

INTERSECTION  ANGLE 

1 

70°  /  50° 

o 

o 

2 

CT> 

O 

O 

O 

o 

20° 

2 

75°  /  50° 

25° 

3 

4^ 

O 

O 

\ 

OJ 

o 

o 

10° 

3 

-p* 

o 

CO 

o 

o 

17° 

3 

50°  /  40° 

10° 

3 

70°  /  47° 

23° 

3 

65°  /  50° 

15° 

Table  5.2:  Ranking  of  best  same  side  stereo  images. 
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5.2.  Distortions  of  Same  Side  Stereo  with  Applied  Squint 
Angles 


The  interpreters  graded  the  stereo  pair  with  50°  and 
70°  elevation  angles  as  best.  Usinq  the  same  input 
parameters  for  the  simulator,  the  squint  anqle  was  varied 
between  0.5°  and  50°.  The  images  obtained  were 
geometrically  restored  to  imaqes  without  squint.  Examples 
are  given  in  Fig.  5.3  (a)-{d). 

The  radiometric  distortions  proved  to  be  "not 
accaptable"  by  applying  a  squint  of  more  than  or  equal  to 
30°.  Again,  the  stereoscopic  impression  is  possible  much 
longer  in  the  northern  and  flatter  part.  The  radiometric 
differences  (illumination,  shadow)  are  much  hiqher  in  the 
southern  part. 

A  summary  of  the  quality  vs.  squint  anqles  is  shown  in 
table  5.3. 

The  front  slopes  were  more  easily  fused  than  the  back 
slopes  where  shadow  effects  add  a  confusion  element  in  the 
viewing  process.  Also  the  limitation  of  the  simulation 
model  which  did  not  take  the  imaqe  resolution  properlv  into 
account,  (described  in  more  detail  in  Appendix  A) ,  may  bias 
the  conclusion  in  favor  of  larqer  squint  anqles  to  be 
acceptable.  Thus  sharp  edges  and  plateaus  were  present  and 
preserved  in  the  simulation. 
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Crossinq  Flight  Tracks  —  Same  Elevation  Angle 


In  the  Greek  Island  SIR-A  picture  imaqinq  with  two 
tracks  crossinq  at  an  angle  of  34°  created  a  qood  stereo 
pair  (grade  "7").  Further  simulations  showed  a  limit  to 
successful  viewing  between  40°  and  50°  of  the  crossinq 
angle.  This  range  of  crossinq  angles  was  investigated  (Fin. 
5.4  and  5.5)  more  closely: 

Successful  stereo  viewinq  was  still  obtained  from 
crossinq  tracks  of  40°  (grade  "5")  ,  but  graded  "1"  for  a 
crossing  angle  of  50°.  The  limit  lies  between  40°  and 
42°,  as  can  be  observed  from  table  5.4.  Trvinq  to  detect 
the  limit  even  more  precisely  was  found  to  he  dependent  on 
the  individual  investigator. 


Figure  5.4:  Flight  and  swatch  configuration 
of  "cross  angle"  stereo. 


CROSSING  UNCLE  50  DEC 


6.  On  the  Accuracy  of  Radar  Stereo  Mapping 


The  theoretical  error  propagation  into  radar  stereo 
model  coordinates  has  been  investigated  in  the  part.  Also 
practical  studies  with  actual  radar  imaqes  have  resulted  in 
an  understanding  of  achievable  accuracies,  The  following  is 
a  short  review  of  recent  work  to  allow  one  to  relate  stereo 
viewability  to  error  propagation  and  accuracy. 


6.1.  Theoretical  Error  Analysis  of  Radar  Stereo  Models 


Error  equations  for  radar  stereo  models  were  derived  bv 
Leberl  (1978,  1979)  for  a  simplified  fliqht  arrangement  as 

shown  in  Figure  6.1.  In  this  local  coordinate  system  the 
vectors  of  sensor  positions  (^'  and  s1')  and  the  velocity 
vectors  of  the  sensor  (s',  s")  are  qiven  as  follows  : 

S'  =  (SX',  Sy',  SZ')  =  (0,  0,  H) 

1'  =  (Sx«,  0,  0) 

8"  =  <Sx’',  sy",  SZ")  =  (0,  B,  A) 

8"  =  (Sx"/  0,  0) 

EO. ’s  6.1 


FLIGHT  LINE 


Figure  6.1.  Coordinate  system 


One  is  dealinq  with  errors  in  sensor  positions  Ms  , 
:1s*  ' )  ,  sensot  velocity  vectors  (ds’,  ds’’)  ami  slant 
rnnqes(dr',dr").  ~ 


(a)  Errors  in  sensor  position  coordinates  : 

It  is  obvious  to  see  that  an  error  dsx’  has  no  effect 
on  y  and  z,  so  that: 

dx(sx')  =  dsx'/2  dx(sx")  =dsx"/2 

<1v(sx’)  =  0  dv(sx*  ' )  =  0 

dz(sx')  =  0  dz(sx'')  =  0 

F,0 .  '  s  6 . 2 


The  effects  of  dsy',  dsv*'  are  less  obvious. 
Figure  6.2  larifies  the  relationship.  There  is  no  effect 
on  x.  In  y,  we  find: 


dy 

1 

with 

1 

dSy'  ' 

So  that  : 

dy  =  H  *  d Sy  '  sin  6' '/( 


H 
r ' 

sin  0 '  • 

sin  (S'  -  9 '  1 ) 

sin  (0'  -  0’  ’)  ) 

EO.  6.3 


But  : 


r '  sin  {90°  -  e • ' ) 
B  sin  (0'  -  0'  ') 


r'  sin  (9'  -  0'')  =  B  sin  {90°  -  0 '  ' ) 


EO.  6.4 


Introducing  EQ.  6.4  in  EQ.  6.3  gives  : 

dy  =  dsy • '  h  tan  0' ’/ B 
so  that,  finally  : 

dy  (sy")  =  d Sy 1  ’  (B  -  v)/B 

EQ.  6 . 6 

We  can  substitute  for  v/B  also  : 

y  y  r'  sin  0'  sin  (90°  +  0'')  1 

B  r'  B  sin  (0'  -  9'')  1  -  cot  0'  tan  0" 
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One  is  dealing  with  errors  in  sensor  positions  (ds', 
ds''),  sensor  velocity  vectors  (ds* ,  ds’’)  and  slant 
ranges  (dr ' ,  dr ' ' ) .  _ 


(a)  Errors  in  sensor  position  coordinates  : 


It  is  obvious  to  see  that  an  error  dsx'  has  no  effect 
on  y  and  z,  so  that: 

dx(sx')  =  dsx'/2  dx(sx '  ')  =  dsx"/2 

dY(sx')  =  0  dy(sx")  =  0 

dz ( sx ' )  =  0  dz ( sx ' ')  =  0 

EO. ' s  6.2 


0" 


* 

I  \ 


Fiqure  6.3.  Errors  due  to  dsz *' -component 
of  sensor  position 


An  error  of  sensor  altitude  does  not  deform  the  model 
y-dimension,  but  merely  shifts  the  model  in  y.  In  heiqht, 
there  is  a  linear  effect  of  model  tilt. 


(b)  Errors  in  sensor  attitude  vectors  : 


Erroneous  sensor  attitudes,  caused  by  errors  of  the 
velocity  vectors  s' ,  s^'  *  r  lead  to  nonintersectinq 
projection  circles.  This  is  illustrated  in  Fiqure  6.4,  for 
the  example  of  an  error  dsv  which  creates  an  antenna 
swinq  k.  We  need  to  define  a  model  point,  that  most 
loqically  is  chosen  halfway  between  the  error-free  point  if 
no  dSy  error  had  been  present. 

From  Fiqure  6.4  we  find  that  an  error  dsv'  affects 
the  x-coordinates  of  the  model  point  P,  shiftinq  it  to 
position  £,  but  leaves  the  v,  z-coordinates  unchanqed,  at 
least  to  a  first-order  approximation  : 
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<3x(Sy«)  *  dSy*  *  y/2  =  dSy*  *  H  *  tan  fi'/2 
dv(Sy')  =  0 

<Jz(sY')  =  0 

dx(sy'*)  —  d  '  *  (y  -B)/2  = 

=  dSy'  '  *  H  *  tan  «"/2 
dy(Sy")  =0 
dz (sv' ' )  =0 

y  EO.'s  6.8 


TOP  VIEW 


Figure  6.4.  Error  due  to  dSy  -  component 
of  velocity  vector 


Similarly,  the  dsz'  dsz ' ' -components  cause  a 
tilt  0  of  the  radar  antenna,  so  that: 


dx(sz ' )  =  dSz'  *  H/2 
dv(sz')  =0 
dz (sz ' )  =  0 

dx(sz")  =  dSz"  *  H/2 
dy(sz")  =  o 
dz ( sz ' •)  =  o 

EO.'s  6.9 


An  erroneous  dsx  '  '  -component  will  not 

affect  the  antenna  attitude  and  thus  not  produce  errors  of 
the  stereo  model.  It  will,  of  course,  affect  antenna 
position  when  integrated  over  time.  This  error  was 
considered  in  Equation  6.2. 
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(c)  Errors  in  slant  ranges: 


A  last  error  of  the  stereo  model  results  from  erroneous 
range  data  r',  r''.  Figure  6.5  shows  that: 

dx(r')  =  0 

dy(r')  =  dr'  *  r'/B  =  dr'(y2  +  H2)1/2/b  = 

=  dr'  *  cos  0''/sin  (0'  -  0'') 
dz(r')  =  dr'  *  r*  *  (y  -  B)/B  *  H  = 

=  dr'  sin  0''/sin  (0'  -  0’’) 

=  dr'  sin  0''/sin  (0'  -  0'') 


dy (r ' ' )  =  -dr1'  *  r''/B  =  -dr''  cos  0'/sin  (0'  -  0'') 

dz  (r  '  ' )  =  -dr''  *  y  *  r”/B  *  H  = 

=  -dr''  *  sin  0'/  sin  (0'  -  0'') 

F.Q. 's  6.10 

These  errors  are  nonlinear  in  the  v  and  z-dimensions. 


^  / 

/ 

/ 

/ 

/ 

/  dy 


Figure  6.5.  Error  due  to  range  error 


6.2  Accuracies  with  Parallax  Measurements 


As  shown  in  Chapter  2,  parallax  difference  measurements 
may  serve  to  compute  heiqht  differences  between  terrain 
points  by  using  the  simplified  equations  (6.11)  for  qround 
range  presentation  or  (6.12)  for  slant  range  presentation : 


dh  =  dpq 

/ 

(  cot 

0'  • 

-  cot 

0') 

EO. 

6.11 

dh  =  dps 

/ 

(  cos 

0  '  ' 

-  cos 

0') 

EO. 

6.12 

The  following  stereo  model  cases  were  used  for  narallax 
measurements: 

(a)  Aircraft  radar  imaqe  pair  of  Granite  Mountain 
(Arizona,  USA)  in  qround  range  presentation. 
This  imaqe  pair  is  shown  in  Fiqure  2.4. 

(b)  Satellite  radar  imaqe  pair  of  Granite  Mountain  in 
slant  range  presentation  with  optical  correlation. 
This  imaqe  pair  is  shown  in  Figure  2.S. 

(c)  Satellite  radar  imaae  pair  of  Los  Anqeles, 

California  in  slant  ranqe  presentation  and  with 
optical  correlation.  This  imaqe  pair  is  shown 
in  Figure  2.6. 

(c)  Satellite  radar  imaqe  pair  of  Los  Anqeles, 

California  in  slant  ranqe  presentation  and  with 
digital  correlation.  This  image  pair  is  shown  in 
Figure  6.6. 


Parallax  differences  between  control  points  as 
differences  of  absolute  parallaxes  were  measured  for  the 
four  radar  stereo  models  usinq  a  stereoscone  and  a  parallax 
bar  and  for  the  two  models  of  Granite  Mountain  also  using  a 
photogrammetr ic  stereoplotter  as  a  comparator.  with  these 
parallax  differences  the  heiqht  differences  between  terrain 
points  were  determined.  The  heiqht  errors  dh  in  control 
points  were  computed  as  the  difference  between  known  heiahts 
h  and  radar-derived  heights  h'. 


B  4 


lei  of  Los  Angeles,  satellite  radar, 
Ltal  correlation 


Since  there  are  systematic  ^  irs  in  the 
radarq rammetr ically  computed  heiqhts  ,  a  correction 
polynomial  was  used  to  eliminate  these  errors: 


h*  = 


3  3 

£  £ 
i=l  3=1 


( i-1)  (i-1) 


EO.  6.13 


One  miqht  expect  that  with  the  use  of  a  flexible 
polynomial  the  systematic  heiqht  errors  can  be  eliminated, 
so  that  only  rerrors  caused  by  measurement  inaccuracies 
remain.  Table  6.1  presents  the  results  of  this  exercise  in 
the  form  of  root  mean  square  residuals  in  the  measured 
control  points.  The  residuals  are  found  in  the  control 
points  between  the  known  heiqht  and  the  polynomial  error 
surface.  The  followinq  conclusions  result: 


The  stereoplotter  was  not  superior  to  simple  parallax 
bar  measurements. 

Systematic  errors  exist  in  all  raw  heiqhts  and  need 
to  be  corrected  with  the  use  of  control  points  and 
correction  polynomials. 


Aircraft  radar  of  Granite  Mountain  provided  hiqher 
accuracies  than  SEASAT  satellite  radar.  This  refers 
to  the  fact,  that  the  aircraft  stereo  model  has  a 
better  stereo  conf iquration,  which  is  expressed  by 
the  relationship  of  stereo  base  and  flvinq  heiaht. 

SEASAT  satellite  radar  of  Los  Anqeles  is  poorer  than 
of  Granite  Mountain  because  of  a  poorer  stereo 
conf iquration  (smaller  base). 

Diqital  and  optical  correlations  led  to  the  same 
performance  fiqures  with  the  imanes  that  were 
employed. 


Imagery  Measurement  Nr.  of  measured  R.m.s. height  error 

points  (m) 
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Table  6.1  Root  mean  square  errors  of  radargrammetric  heights  after 
polynomial  correction 


6.3  Riqorous  Radar  Stereo  Mapping 


6.3.1  Mathematical  Formulation 


For 
for  each 


each  qround  point  there  are  two  condition  equations 
imaqe : 


(1) 

The 

ranqe 

cond i tion : 

1  £ 

-  s  | 

=  r 

(2) 

The 

squint 

anqle  condition: 

sin 

t  =  s 

7  (p  -  s)  /_(  s  j 

EO.  6.14 


eo.  tt.it* 


For  the  formation  of  a  radar  stereo  model  we  need  to 
use  a  pair  of  these  conditions: 


r ' 

r  "  - 


0 
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s-  *  (  p 

sint'  s'_|  * 

s"*  (  p 
sin  f ' ’ *  7  s ’ '  I  * 


s'J 

p  -  s'  |_=  0 
s' 7)  - 

p  -  s  ’  ’  I  =  0 


EO.  6.16 


are 


The  sensor  positions  s  and  the  sensor  veloc 
represented  as  polynomials  of  the  imaqinq 


itv  vectors 
times  t: 


s  = 

I  = 

(SX  t  Sv  r  Sz  )  -  S 
(SX  t  Sx  i  Sx  )  =  S 

% 

— 

EO.  6.17 

Imaq ing 

times  t  can  be 

expressed 

bv 

radar  imaqe 

coordinate  x 
as  follows: 

and  slant  range  r 

by  radar 

imaqe 

coordinate  v 
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X 

*  x  +  cx 

i-s 

ii 

3 

< 

*  v  +  cv 

EO.  6.18 
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The  x-axis  of  the  radar  image  system  is  pointing  in 
flight  direction  of  the  imaging  sensor.  mx  and  mv  are 
image  scale  factors,  cx  is  the  time  corresponding  to  image 
coordinate  x  =  0  and  Cy  is  the  sweep  delav. 

If  we  assume  error-free  sensor  orientation  and  image 
parameters  (mx,  my,  cx,  Cy) ,  the  equations  6.16 
contain  the  unknowns  X,  Y,  Z  of  the  object  point  and  the 
measurements  x',  y',  x",  v"  for  the  image  coordinates  of  the 
left  and  the  right  radar  image.  Linearization  leads  to: 


C  *  v  +  D  *  dp  +  w  =  0 


EO.  6.19 


C  and  D  are  coefficient  matrices  of  the  linearization, 
the  vector  v  contains  the  corrections  for  the  image 
coordinate  measurements,  the  vector  dp  contains  increments 
for  approximately  known  ground  coordinates  of  the  object 
point  and  the  elements  of  vector  w  are  the  contradictions  of 
Equ.  6.16  when  introducing  image  coordinate  measurements 
and  approximate  X-,  Y-,  Z-  values.  The  unknowns  do  can  be 
found  by  a  least  squares  adjustment  with 


djD  =  -  (DT  (  C  G  cV1  dt  (  £  G  C  )  1  W 

EO.  6.20 

If  orientation  and  image  parameters  are  not  error-free, 
they  may  be  improved  by  a  radar  block  adjustment  or  one  can 
use  ground  control  points  to  determine  calibration 
polynomials. 

In  the  case  of  calibration  the  X,  Y,  Z  -  coordinates  of 
the  ground  control  points  are  converted  to  left  and  right 
radar  image  coordinates  (xt',  vt',  xt",  yt")  usinq  the 
erroneous  sensor  parameters.  The  difference  between  these 
transformed  and  the  measured  image  coordinates  is 
approximated  with  correction  polynomials: 
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so  that 

xt'  =  x'  +  dx' (x' ,  y') 
yt'  =  y'  +  dy' (x* ,  y' ) 
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6.3.2  Results  Obtained  with  Rigorous  Stereo  Radararammetrv 


The  procedure  described  above  to  convert  radar  imaqe 
coordinates  x',  y',  x",  v"  to  object  coordinates  was  applied 
to  a  SEASAT  radar  stereo  model  of  the  area  of  Los  Anqeles 
and  to  a  SIR-A  radar  stereo  model  of  the  qreek  islands 
Cephalonia  and  Ithaka.  The  SEASAT  radar  stereo  model  is 
shown  in  Figure  6.7  and  the  SIR-A  model  previouslv  was 
presented  in  Chapter  2,  Figure  2.13. 


Image  coordinates  were  measured  with  the  analytical 
plotter  KERN  DSR-1  used  as  stereo  comparator  and  transformed 
to  object  coordinates  Xtf  yt,  zt.  These  computed 
ground  point  coordinates  then  were  compared  with  the  given 
values  X,  Y,  Z.  Table  6.2  summarizes  the  results  that  were 
obtained  with  the  mentioned  algorithm.  Different  kinds  of 
correction  polynomials  were  used  to  calibrate  the  erroneous 
sensor  data. 

Since  the  radar  image  parameters  of  these  models  are 
unknown,  at  least  a  calibration  polynomial  of  linear  order 
in  x  and  y  is  neccessarv  to  qet  reasonable  results.  while 
the  height  accuracy  of  the  SEASAT  stereo  model  comoares  well 
with  that  of  SIR-A,  it's  accuracy  in  olanimetrv  is  somewhat 
poorer.  This  refers  to  the  steeper  look  angles  of  the 
SEASAT  imaging  system  (approximately  22  deqrees) ,  where  for 
example  an  error  in  slant  range  causes  larqer  errors  for 
planimetry  than  for  heights  (note  from  chanter  7.1  that  the 
planimetric  error  is  affected  bv  the  cosine  and  the  height 
error  by  the  sine  of  the  elevation  anqle) .  Since  the  look 
angles  for  the  SIR-A  radar  system  are  about  45  deqrees,  also 
the  errors  in  planimetry  are  of  the  same  order  of  m^qnitude 
than  those  for  heights. 

Further  we  see  that  with  the  use  of  flexible 
calibration  polynomials  of  second  order  in  x  and  v  the 
results  may  be  improved.  However,  the  improvement  for 
SEASAT  stereo  is  about  30  percent  while  for  SIR-A  stereo  it 
is  only  about  20  percent.  This  is  not  as  much  as  one  miqht 
have  expected. 

Also  the  use  of  a  poorer  but  well  distributed  qround 
control  density  leads  to  quite  reasonable  results  that  are 
not  much  poorer  than  those  obtained  with  high  qround  control 
density. 


i  g  u  r  e  b  .  7  SEA5AT  stereo  model  oi  Los  Angeles,  California, 
same  side  geometry 
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7.  CONCLUSIONS  AND  OUTLOOK 


The  study  concerns  stereoscopic  viewinq  of  side- look inq 
radar  images.  Current  understanding  of  this  effect  is  based 
on  theoretical  considerations  of  radar  parallax,  of  error 
propagation  and  work  with  actual  radar  imaqes.  This, 
however,  is  too  limited  to  be  satisfactory:  the  precise 

capabilities  and  limitations  of  radar  stereoscopy  remained 
unclear . 

In  the  framework  of  the  current  study  one  therefore 
developeded  a  technique  for  computer  simulation  of  radar 
images.  Using  this  one  can  vary  a  number  of  significant 
parameters  of  the  radar  imaging  process  so  that  a  better 
understanding  is  obtained  of  radar  stereoscopic  effects. 

This  report  presents  the  startinq  point  for  the  studv 
by  a  review  of  basic  principles  of  radar  stereo- imaqing, 
parallax  and  vertical  exaggeration.  It  then  proceeds  to 
describe  the  chosen  technique  of  radar  image  simulation. 
Emphasis  is  on  radar  image  geometry  since  stereoscopy  is 
largely  a  geometric  effect.  An  image  space  algorithm  was 
chosen  as  an  appropriate  method  of  simulation  comou tat ions. 

The  validity  of  computer  generation  of  imaqes  was  with 
simple  geometric  shapes  and  with  various  real  radar  imaqes 
of  areas  where  digital  terrain  elevation  models  exist.  A 
great  number  of  images  was  then  generated  with  various  look 
angles,  squint  angles  and  fliqht  directions  to  allow  for  a 
complete  review  of  stereo  viewinq  limitations.  'J’he 
conclusions  are  as  follows: 

(a)  Same-side  stereo  from  parallel  flight  lines  produces 
the  most  distinct  vertical  exaqqeration  at  steeoest 
look  angles,  e.g.  10°  and  20°  off-nadir,  with 
intersection  angles  of  only  10°. 

(b)  This  same  side  geometry  allows  successful  stereo- 
-viewing  with  intersection  angles  of  up  to  60°, 
namely  look  angles  off-nadir  of  20o  and  of 

80o,  in  the  extreme  case. 

(c)  Larger  differences  of  look  angles  cannot  be  fused 
stereoscopically  due  to  excessive  shadowinq  in  one 
image  combined  with  excessive  lay-over  in  the  other 
image. 

(d)  Opposite-side  stereo  from  parallel  flight  lines  was 
not  usable  for  stereo  viewing.  However,  unprocessed 
raw  images  with  oriqinal  qrav  values  were  used  for 
analysis. 
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(e)  Crossing  flight  line  stereoscopy  is  feasible  up  to 
intersection  angles  of  about  42°. 

(f)  The  use  of  squinted  radar  images  does  not  destroy 
the  stereoscopic  effect  if  squint  anqles  amount  to 
a  total  of  30o  or  less. 

As  a  by-product  of  the  study  also  a  mathematical  model 
was  developed  for  parallax-computations  with  crossinq 
flight-line  stereo.  An  understanding  was  obtained  For  the 
occurrence  of  x  and  y  -  parallaxes  and  for  the  required 
orientation  of  the  eye  basis.  An  analysis  produced  maximum 
y  parallaxes  for  one  existing  340  crossinq  anqle  SIR-A 
stereo  pair  in  the  amount  of  160  m  on  the  ground  (SO  km 
swath) . 

Also  a  number  of  rigorous  radarqr ammetr ic  computations 
were  performed  on  actual  radar  image  pairs  from  both 
aircraft  and  satellite,  with  accuracies  ranginq  from  4R  m 
(aircraft)  to  143  m  with  Seasat  satellite  data.  SIR-A 
cross-angle  stereo  resulted  in  78  m  to  123  m  heiqht  errors, 
depending  on  the  ground  control  used. 

The  resources  for  the  study  did  not  allow  to  cover  all 
aspects  of  the  area  of  radar  stereoscopy,  nor  of  the  area  of 
topographic  shape  reconstruction  from  radar  imaqes. 

A  next  and  important  step  is  the  addition  of  both 

(a)  ground  resolution, 

(b)  thematic  variation 

in  the  computer  simulation.  The  variability  of  ground 
resolution  will  not  affect  the  "viewabilitv"  of  radar  stereo 
models,  but  it  will  affect  the  surface  definition  and 
therefore  the  accuracy  of  measurements.  This  is  a 
significant  element  in  the  evaluation  of  radar  stereo  work. 

Thematic  variability  leads  to  the  incorporation  of 
texture,  use  of  several  backscatter  curves  and  a  chance  to 
more  realistically  simulate  the  radar  imaqes.  It  is  also 
meaningful  to  better  understand  the  relative  interdependence 
of  height  and  planimetric  measurement  accuracies. 

Clearly  the  work  would  benefit  from  studvinq  a  wide 
range  of  terrain  types.  This  was  not  possible  in  the 
current  effort  due  to  limitations  in  both  time  and 
resources. 

The  improved  capability  of  radar  imaqe  simulation 
should  be  combined  with  current  radar  imaqe  stereo 
measurement  techniques  on  photoqrammetr ic  analytical 
plotters.  This  will  then  close  a  complete  loop  to  obtain  a 


topoqraphic  stereo  mappinq. 


Appendix  A 


Refined  Description  of  SIMRISA 

Fig.  A. 1  shows  a  table  of  contents  of  the  proqram 
package.  Each  rectangle  defines  a  program  unit;  program 
functions  like  input  routines  or  help  functions  do  not 
appear  in  the  table  of  contents,  but  are  added  after  the 
description  of  the  main  routines. 

Documentation  time  during  and  throughout  the 
development  cycle  was  unfortunately  -  like  always  -  rather 
limited.  But  as  good  documentation  is  essential  to  the 
continued  success  of  a  program  development  effort,  our 
choice  was  a  HIPO-like  (Hierarchy  plus  Input-Process-Output) 
documentation  design.  This  has  the  advantage  of  not 
describing  the  program  logic  -  or  at  least  not  until  a  very 
low  level  -  but  the  functions  of  the  system.  Thus  changing 
programming  logic  during  the  implementation  did  not  result 
in  a  redesign  of  a  whole  diaqram  or  flow  chart. 

The  next  pages  show  a  short  documentation  of  the  main 
routines.  Each  routine  expects  some  input  and  creates  some 
output,  the  corresponding  process  is  described  with  the  help 
of  programming  statements,  normal  English  words  and 
structograms.  It  should  be  easy  enough  to  read. 
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